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Abstract 

The relativistic transformation properties of the heavy quark-antiquark inter- 
action potential are considered in the framework of the relativistic quark model. 
A special attention is paid to the long-range (confining) contribution to the spin- 
independent part of qq interaction. The retardation effects are consistently taken 
into account. 

The relativistic properties of the quark-antiquark interaction potential play an impor- 
tant role in analysing different static and dynamical characteristics of heavy mesons. The 
Lorentz-structure of the confining quark-antiquark interaction is of particular interest. In 
the literature there is no consent on this item. For a long time the scalar confining kernel 
has been considered to be the most appropriate one [0], ||. The main argument in favour 
of this choice is based on the nature of the heavy quark spin-orbit potential. The scalar 
potential gives a vanishing long-range magnetic contribution, which is in agreement with 
the flux tube picture of quark confinement of ||, and allows to get the fine structure 
for heavy quarkonia in accord with experimental data. However, the calculations of elec- 
troweak decay rates of heavy mesons with a scalar confining potential alone yield results 
which are in worse agreement with data than for a vector potential j|, [5]]. The radia- 
tive Ml-transitions in quarkonia such as e. g. J/ip —>■ r} c ^f are the most sensitive for the 
Lorentz-structure of the confining potential. The relativistic corrections for these decays 
arising from vector and scalar potentials have different signs j|, || . In particular, as it has 
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been shown in ref. agreement with experiments for these decays can be achieved only 
for a mixture of vector and scalar potentials. In this context, it is worth remarking, that 
the recent study of the qq interaction in the Wilson loop approach || indicates that it can- 
not be considered as simply a scalar. Moreover, the found structure of spin-independent 
relativistic corrections is not compatible with a scalar potential. A similar conclusion has 
been obtained in ref. J7| on the basis of a Foldy-Wouthuysen reduction of the full Coulomb 
gauge Hamiltonian of QCD. There, the Lorentz-structure of confinement has been found 
to be of vector nature. The scalar character of spin splittings in heavy quarkonia in this 
approach is dynamically generated through the interaction with collective gluonic degrees 
of freedom. Thus we see that the spin-dependent structure of (qq) interaction is well 
established now, while the spin-independent part is still controversial in the literature. 

In preceding papers ||, |§ we have developed the relativistic quark model with the (qq) 
potential consisting of the perturbative one-gluon exchange part and a nonperturbative 
one which is a mixture of the Lorentz scalar and vector confining potentials: 



where k = p — q, D^ v is the gluon propagator in the Coulomb gauge and T M is the effective 
vector long-range vertex, containing both the Dirac and Pauli terms 



we can treat the parameter (1 + k) as the nonperturbative (long-range) chromo magnetic 
moment of the quark and k as its anomalous part (flavour independent). 




(1) 



u a ,b(p) are the Dirac bispinors. Using the identity 



(2) 





where 




In order to reproduce the linear confining potential 



Konf(r) = Ar + B 



in this limit one should put 





(l-e)(Ar + B) 



(4) 
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where e is the mixing parameter. 

Now assuming that both quarks are heavy enough we evaluate the (v 2 /c 2 ) relativistic 
corrections to the static potential ([3D, (f|). For the one-gluon exchange part the retardation 
effect is taken into account by the contribution of transverse gluon exchange and thus it 
is sufficient in the adopted approximation to set k° = 0. For the confining part one 
should utilize a different procedure (see |TI], [II], The Fourier transform of the linear 
potential Ar in the momentum space looks like: 

A J d 3 rre" ikr = ~ A KTi, k = p - q. (5) 

The natural (though not unique) relativistic extension (dependent only on the four- 
momentum transfer) of expression (|) is to substitute (— k 2 ) — > (k^ — k 2 ) and thus 



f f 



(fc 2 -k 2 ) 2 ' 



(6) 



Now as mentioned above we should choose the procedure of fixing ho- On the mass shell 
due to energy conservation we have k Q = 0. So ko may be considered as the measure 
of deviation either from the mass shell or from the energy shell. We choose the second 
possibility and set ko equal to e a (p) — e a (q) or to e&(q) — e&(p). Then in the symmetrized 
form |T| we get 

fc 2 = -( ea (p)-e a (q))(e 6 (p)-6 fc (q)), (7) 



e»MP) = VP + m a 



This form is not unique and other possible expressions for k$ are discussed in [ |I2] , |TB| ] . In 



favour of choice (0) we mention the following arguments. It is well-known Jl0|, [14j] that for 
the one-photon exchange contribution in QED only choice (|7]) in the Feynman (diagonal) 
gauge leads to the same correct result (the Breit-Fermi Hamiltonian) as the prescription 
ko = in the Coulomb (or transverse Landau) gauge. The same is naturally true for the 
one-gluon exchange contribution in QCD. Moreover as shown in ref. [Ill] for any effective 



vector potential generated by a vector exchange and its couplings to conserved vector 
currents (vertices) there is the so-called instantaneous gauge which plays the role of the 
Coulomb gauge. In the instantaneous gauge the prescription k = reproduces the same 
result as the expansion in k\ fixed by eq. (|7|) in the diagonal gauge used here. The other 
reason to utilize prescription (0) is the reproduction of the correct Dirac limit in this case 

g. 

The (p 2 /m 2 ) expansion of (Q) with the account of eq. (0) yields: 

\ _ (p 2 - q 2 ) 2 



2m a mJk| 2 



[Mp) - £«(q))Mp) - e»(q)) + k 2 ] 2 |k 

= ^-2d^ l(k ' P)2 + 2(k ' P)(k ' q) + (k ' q)21 ' <8) 
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After the Fourier transform of eq. (§D we obtain in the configuration space: 



-8nA 



d*k ^. r J_ 



(2tt 



1 - 



(P 2 ~ q 2 ) 2 
2m a m b \h\ 2 



At - 



Ar 



2m a m b 



p 2 + ^ 



(9) 



where the notation {. . .}w means the Weyl ordering prescription for p and r variables. 

Now we turn to the constant B term in eq. (||). The Fourier transform of it gives in 
the momentum space 



B J d 3 re" ik - r = B(27i) 3 5(k), k = p - q. 



(10) 



The simplest relativistic version of eq. (10) is multiplying it by the energy factor e(p)/m, 
which in the symmetric form looks like 



m a m b 



<5(p-q). 



Expanding eq. (|TTD in (p 2 /m 2 ) we get 



mi 



(12) 



(13) 



All other relativistic corrections of order (p 2 /m 2 ) have been considered in ref. ||. After 
using the Weyl ordering notations, expressions (9) of ref. || take on the form: 



<J(p - q) 

So in the configuration space the constant term acquires the form 

B 
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\ 




m 2 b ) 



V (r) + 
1 

-i 

2m a m b 



1 f 1 



m 2 a + m\ t 



V ( 



Coul 



2 IP ' r 



A^coui(r) + 

2 



1 ( l+2K a | 1 + 2K b 



mt 



mi 



^""a""6 ^ L J ) W 



Adding contributions ([|) and (p~3|) to the above expressions we obtain the complete spin- 
independent part of the (qq) potential (n a = n b = k): 



Vsi(r) = V (r) + V VD (r) + 



1 ( 1 



m 2 ^ ml 



A[Vboui(r) + (l + 2«)Kr»], (14) 



where Vo(r) is given by eqs. 
Vvn(r) ' 



2m a m b 



4 a 
3 r 



. For the velocity-dependent part Vyr>(r) we have 

2" 



P 2 + ^ 



W 
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m a m b 



(15) 



Now representing eq. in the form 

V YD (r) 1 



pV 6c (r) + 



m a m fc 
A + |p 2 ^e(r) 



-K(r) 



p-r 



mi 



-K(r) 



(16) 
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with 



V bc (r) 
V de (r) 
VJr) -- 



-£ + (*-,)* + (l-e)2I; 

- £ -Ar+( 1 -- £ -)B; 
2 U 2/ 

- — - -Ar; VJr) = 0, 
3r 2 ' v ' 



(17) 



we are able to test the fulfilment of the exact Barchielli, Brambilla, Prosperi (BBP) rela- 



tions [15], which follow from the Lorentz invariance of the Wilson loop. In our notations 
these relations look like 



de 



-V bc + -^ = 
2 4 



(in the original version V bc = —V b — \ V C and Vd e = Vd + \ V e ). One can easily find that the 
functions (|i~7D identically satisfy relations ( |i~8|) independently of values of the parameters 
e and k. This is a highly nontrivial result. For the perturbative one-gluon-exchange part 
of VVd our expressions for V b , . . ., V e are the same as in [[l|, |i~6f , but for the confining 
(long-range) part they are different, namely the result of refs. [fT5" , |16| (from the minimal 
area law) is as follows: 



V bc (r) 
V d Jr) 



2a 



1 



+ -Ar; 
3r 6 

— Ar - -B- 

6 4 ' 



V c (r) 
V, = 



_2a ± 

3r 

--Ar. 
6 



1 , 

-Ar: 
6 ' 



(19) 



No value of e in eqs. (|i~7|) can reproduce the above result. The terms with the Laplacian in 
(|I4]) coincide only for k = and e — 0, i. e. for purely vector confining interaction without 
the Pauli term in the vertex (|2|). Our expressions ( |14|) and (|T5"|) for purely vector (e = 0) 
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coincide with those of ref. 13 



except 

(scalar potential) is the same as in fT5| . 

-B/2 



and purely scalar (e = 1) interactions and for k 
for the constant B term. Our B term for e = 1 

The B term from ref. [13| does not satisfy the BBP relations (it gives contribution 
only to Vde)- Our result for the scalar (e = 1) confining potential also differs from 
the one obtained in ref. fTTH , where the prescription k = was used and as a result the 
contribution of retardation was lost. The differences between our results and the results 
presented in ref. |18| originate from the use of specific models such as minimal area law, 
flux tube, dual superconductivity and stochastic vacuum. 

The spin-dependent part of the (qq) potential is given in ref. (n a = Kb = k)' 



^SD(r) 



-T-Vboui(r)L ■ (s + s b ) 



m a nib r dr 
1 1 d 
2 m 2 r ^ r 

1 1 d 



{[Wr) - Konf(r)] + 2(1 + K) (l + ^) Konf( r )} L ■ Sa 

[Vboui(r) - Konf(r)] + 2(1 + K)(l + —) 



2m\ r 

<l2 [Vcoui + (1 + «) 2 C nf (r)] 



dr 



^(s a • r)(s a • r) - s a • s 6 



A[V C oui(r) + (1 + K)Vl(r)]s B ■ s b 



3m a m b 



(20) 



where L is the orbital momentum and s a ^ are the spin momenta. 

The correct description of the fine structure of the heavy quarkonium mass spectrum 
requires vanishing of the vector confinement contribution. This can be achieved by putting 
1 + k — 0, i.e. the total long-range quark chromomagnetic moment equals to zero, which 
is in accord with the flux tube || and minimal area 18 1 models. One can see from 
eq. (|2"0D that for the spin-dependent part of the potential this conjecture is equivalet to the 
assumption about the scalar structure of confinement interaction |]]. The specific value 
of vector-scalar mixing parameter e — — 1 provides the correct description of radiative 
decays of heavy quarkonia . 

In this way setting k — —1, we obtain: 



1 /4a 
2 



3r 3 
1 



.4 

r 

4a 



m a mb V3r 3 
1 /4a s 



m a mb \3r 3 



L • (s a + s b ) 



S7T 



/4a s 



3m a m b V 3 

3 

— (s a ■ r)(s& ■ r) - s a ■ s b 



S(r)s a ■ s b 



(21) 



re 



The 



Expression fl2~T]) for Vsd completely coincides with the one found in refs. 
Gromes relation is identically fulfilled. Our result supports the conjecture that the long- 
range confining forces are dominated by chromoelectric interaction and that the chro- 
momagnetic interaction vanishes. It is also in accord with the dual superconductivity 
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picture |19], yj|. It is important to mention, that our relativistic quark model is in com- 
plete agreement with the heavy quark effective theory (HQET). The model correctly 
reproduces HQET results for heavy-to-heavy and heavy-to-light weak transition matrix 
elements with the specific choice of the parameters e = — 1, k = — 1 (see ref. |]] for details) 
in accord with the ones found previously p0| . It includes the proper description of 



invariant form factors within the inverse heavy-quark-mass expansion up to the terms 
of order of 1/rriQ. The mass spectra of D and B mesons have been also calculated in 
our model in complete agreement with the HQET predictions and available experimental 



data It is interesting to note that the relations which are equivalent to the BBP re- 
lations fll8D can be obtained by use of the reparametrization invariance (in four-velocity) 
within HQET [E2|. The phenomenological implications of retardation corrections will be 
considered elsewhere. 
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